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A time reversal invariant (TRI) topological superconductor has a full pairing gap in the bulk and 
topologically protected gapless states on the surface or at the edge. In this paper, we show that in 
the weak pairing limit, the topological quantum number of a TRI superconductor can be completely 
determined by the Fermi surface properties, and is independent of the electronic structure away 
from the Fermi surface. In three dimensions (3D), the integer topological quantum number in a 
TRI superconductor is determined by the sign of the pairing order parameter and the first Chern 
number of the Berry phase gauge field on the Fermi surfaces. In two (2D) and one (ID) dimension, 
the Z2 topological quantum number of a TRI superconductor is determined simply by the sign of 
the pairing order parameter on the Fermi surfaces. We also obtain a generic and explicit expression 
of the Z2 topological invariant in ID and 2D. 
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Since the discovery of quantum spin Hall effectii^i^ii^, 
topological insulators (TI) in both two dimensions (2D) 
and three dimensions (3D) have generated great inter- 
est both theoretically and experimentally^'^iii^'^. Fol- 
lowing this development, recent attention has focused 
on time reversal invariant (TRI) topological supercon- 
ductors and superfluid&iSiiiii^ii^. There is a direct anal- 
ogy between superconductors and insulators because the 
Bogoliubov-de-Gennes (BdG) Hamiltonian for the quasi- 
particles of a superconductor is analogous to the Hamilto- 
nian of a band insulator, with the superconducting gap 
corresponding to the band gap of the insulator. "^He- 
B is an example of such a topological superfluid state. 
This TRI state has a full pairing gap in the bulk, and 
gapless surface states consisting of a single Majorana 
coneiSiiiiiSJ^. In fact, the BdG Hamiltonian for ^He-B 
is identical to the model Hamiltonian of a 3D topological 
insulator proposed by Zhang et a£. In 2D, the classifica- 
tion of topological superconductors is very similar to that 
of topological insulators. Time-reversal breaking (TRB) 
superconductors arc classified by an integer— li^, similar 
to quantum Hall insulators^^, while TRI superconduc- 
tors are classifie d^°i^^'^^i^^ by a Z2 invariant in ID and 
2D, but by an integer (Z) class in 3D. The integer-valued 
topological invariant in 3D can be written as a winding 
number over the entire momentum spacei^, while an ex- 
plicit expression of the Z2 topological invariants in ID 
and 2D has not been obtained in the literature. 

Despite the similarity of topological insulators and 
topological superconductors, there is a key physical dif- 
ference between them. Starting from a Fermi liquid "nor- 
mal state" with some attractive interaction, supercon- 
ductivity is induced at low enough temperatures due to 
the Cooper instability of the Fermi surface. At least right 
below the transition temperature. Cooper pairing is only 



important around the Fermi surface, so the topological 
properties of such a superconductor are determined com- 
pletely by the properties in the neighborhood of the Fermi 
surface, rather than that of the full Brillouin zone, as in 
the case of a topological insulator. 

Motivated by such an observation, in this letter we 
obtain simple and explicit physical criteria for TRI topo- 
logical superconductivity in 1, 2 and 3 dimensions in the 
weak pairing limit, where the pairing is only important 
in a small neighborhood of the Fermi surface. In 3D, the 
Fermi surface topological invariant (FSTI) of a TRI su- 
perconductor is determined by the sign of the pairing or- 
der parameter and the first Chern number of each Fermi 
surface. Here the first Chern number of a Fermi surface 
is defined by the net flux of the Berry phase gauge field 
penetrating the Fermi surface. This is quantized as long 
as the Fermi surface is a smooth two-dimensional man- 
ifold. Based on this Fermi surface criterion for the 3D 
topological superconductor, the Z2 FSTI's in ID and 2D 
can be obtained by dimensional reductioni^. The crite- 
ria for a Z2 nontrivial superconductor in the weak pairing 
limit is simple: a TRI superconductor is nontrivial (triv- 
ial) if there are an odd (even) number of Fermi surfaces 
with a negative pairing order parameter. For example, 
the superconductivity in a 2D Rashba system is nontriv- 
ial if the pairing on the two Fermi surfaces has opposite 
sign. Inspired by such Fermi surface formulas, we also 
obtain an explicit expression of the Z2 FSTI in ID and 
2D which applies to generic superconductors beyond the 
weak pairing limit. 

Topological invariant in 3D TRI superconduc- 
tors. We start from a generic mean-field Hamiltonian 
of a 3D TRI superconductor, which can be written in 
momentum space as^^ 
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In general, V'k is a vector with components and ft-k and 
Ak are Nx N matrices. The matrix T is the time-reversal 
matrix satisfying T^hk.T = h^^, T'^ = -I and T^T = I. 
We have chosen a special basis in which the BdG Hamil- 
tonian has a special off-diagonal form. It should be 
noted that such a choice is only possible when the sys- 
tem has both time-reversal symmetry and particle-hole 
symmetry. These two symmetries also require TA]^ to 
be Hermitian, which makes the matrix /ik + i'^Aj^ gener- 
ically non-Hermitian. The matrix /ik + iTI\\^ can be 
decomposed by singular value decomposition (SVD) as 
/ik + iT/^l^ = C/k-OkVk with [/k, Vk unitary matrices and 
_Dk a diagonal matrix with non-negative elements. One 
can see straightforwardly that the diagonal elements of 
-Dk are actually the positive eigenvalues of H\^. For a 
fully gapped superconductor, is positive definite, and 
we can adiabatically deform it to the identity matrix I 
without closing the superconducting gap. During this 
deformation the matrix /ik + *^Aj^ is deformed to a uni- 
tary matrix Qk = ^k^k G U(N). As shown in Ref.^^, the 
integer-valued topological invariant characterizing topo- 
logical superconductors is defined as the winding number 
of Qk: 



with e*^"k = (e„k + *<5„k) / krik + i'^nkl- In tbc weak pair- 
ing limit, we take i5„k to be nonzero only in a small neigh- 
borhood —e<E<€oi the Fermi level. As shown in Fig. 
[U the phase 0„k changes from to ±7r across the Fermi 
level, with the sign determined by the sign of (5„k- In the 
limit e — > 0, such a domain wall configuration of 0„k can 
be expressed by the formula 



-7rv„kSgn ((5„k) 5 (e„k) 



(5) 



in which v„k — VkCnk is the Fermi velocity. It should 
be noted that for a gapped superconductor (5„k remains 
nonzero for all k on the Fermi surfaces, so the sign of (5„k 
is fixed on each Fermi surface. 
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Now we study Qk in the weak pairing limit. For sim- 
plicity, in the following, we will assume the Fermi surfaces 
are all non-degenerate, and there arc no lower dimen- 
sional zero-energy defects such as point or line nodes. All 
our conclusions can be easily generalized to more generic 
cases. When the Fermi surfaces are non-degenerate, and 
the weak pairing term Ak is only turned on around the 
Fermi surfaces, the matrix elements of T\\^ between dif- 
ferent bands are negligible. Thus, to leading order we 
have 

/ik + i'?'Aj^ ~ ^(e„k + i(5„k) |«, k) (n,k| (3) 

n 

with(5„fc EE (7i,k|TAj^ |n,k) e M 

where |n,k) are the eigenvectors of /ik. Physically, (5„k 
is the matrix element of A]^ between |n, k) and its time- 
reversed partner |n, — k) — |7i,k). In this approxima- 
tion, the matrix Qk is given by 



Qk 
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FIG. 1: (a) The path of e„k -I- iSnk in the complex plane 
for positive (red) and negative (blue) (5„k around the Fermi 
surface, (b) ^„k and tnk vs momentum k. The change of 
S„k across kp is — tt (-I-tt) when (5,ik is positive (negative), as 
shown by the red (blue) curve. 



Once the behavior of Onki in the Brillouin zone is sim- 
plified to Eq. ([5]) in the weak pairing limit, the wind- 
ing number ^ can be simplified to the following simple 
FSTI: 



Nw = 



i^sgn((5s)Cis 



(6) 



where s is summed over all disconnected Fermi surfaces, 
and sgn(5s) denotes the sign of (5„k on the s-th Fermi 
surface. Cis is the first Chern number of the s-th Fermi 
surface (denoted by FSs): 



Cis^^ I dn^' (9^a«,(k) - 5,a..(k)) 



(7) 



with asi = —i{sk.\d/dki\s\i) the adiabatic connection 
defined for the band |sk) which crosses the Fermi surface. 
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and drJ*^ the surface element 2-form of the Fermi surface. 
More details of the derivation of Eq. ^ arc included in 
the supplementary materiaU^. 

As an example, consider a two-band Hamiltonian /ik — 
/2m — n + ok ■ a For /i > 0, the system has two Fermi 
surfaces which are concentric spheres around k = 0. 
(The two-band model should be regularized on the lat- 
tice, but the lattice regularization is unimportant as long 
as no other Fermi surfaces are introduced.) Denoting 
the electron states at the inner (outer) Fermi surface by 
|k, -l-(— )), we have cr • k |k, ±) = ± |k| |k, ±). It is easy to 
check that the two Fermi surfaces carry opposite Chern 
number C± = ±1. Thus, according to we can obtain 
a topological superconductor if the two Fermi surfaces 
have opposite signs of pairing. The time-reversal matrix 
is T = iuy in this system. If we have Ak — lAoffy, then 
zTAjjj. = AqI which has the same sign on the two Fermi 
surfaces and leads to Nw = 0. On the other hand, if we 
have Ak = iAoaya • k, then iTAj^. = Aqct • k has oppo- 
site sign on the two Fermi surfaces, so that Nw = 1 if 
Aq > 0. If we take the limit a ^ 0, we obtain a topolog- 
ical superconductor with quadratic kinetic energy term 
and pairing Ak = iAoaya ■ k, which is exactly the BdG 
Hamiltonian of the ■^HeB phase. This example also illus- 
trates how the FSTI ([6]) can be generalized to systems 
with degeneracies on the Fermi surface: One can always 
add a small perturbation proportional to TA]^ to the 
Hamiltonian to lift the degeneracy, while preserving the 
topological properties of the superconductor. 

Dimensional reduction to 2D. The FSTI can be 
generalized to lower dimensions i.e. 2D and ID. The TRI 
topological superconductors in 2D and ID are related to 
the one in 3D bv dimensional reduction, similar to the 
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FIG. 2: Dimensional reduction from a 3D TRI superconduc- 
tor to a 2D TRI superconductor. The 2D TRI superconductor 
corresponds to the 6 — tt section of a 3D superconductor. The 
Fermi surfaces with blue (red) color are those with positive 
(negative) pairing amplitude 5s. 

Due to the same symmetry reason as the 3D case, 
the BdG Hamiltonian ffk of a 2D TRI superconductor 
can also be written in the form of ([T]), so that one can 
also define a matrix Qk G U(N) for the 2D case. Since 
H2([/(iV)) = 0, we can always find a smooth deformation 
Qk,o, S G [0; ■"'] which interpolates between Qk and the 



identity I: 

Qm = {q;^Jz° (8) 

It should be noted that Qk satisfies T^QitX = Q^k '^^^ 
to time- reversal symmetry. Thus, if we define Qk,-e — 
T^Q^i^ gT for 6 e [0,7r], we obtain Qk,^ for 6 e [-tt, tt] 
which is continuous and periodic in 9 -\- 2it. Consid- 
ering as a momentum in an additional dimension, Qk^g 
describes a 3D TRI superconductor, which is character- 
ized by the winding number ([2]). If there are two differ- 
ent interpolations Qk.e and g which both interpolate 
between Qk and I, it can be shown that time-reversal 
symmetry requires their winding numbers to be differ- 
ent by an even number: Nw{Q) — Nw{Q') = mod 2. 
Thus the parity (~1)^" (Q) is independent of the choice 
of interpolation path, and is a Z2 topological invariant 
uniquely determined by Qk- 

Now we study the expression of such a Z2 invariant 
in the weak pairing limit. In this limit, the interpola- 
tion of Qk to Qk,e is equivalent to interpolating the ID 
Fermi circles of the 2D normal state Hamiltonian /ik to 
Fermi surfaces in a 3D Brillouin zone parameterized by 
{kx,ky, 9). We can simply extrapolate the pairing on the 
Fermi circles to the Fermi surfaces, as illustrated in Fig. 
([2]). If the Fermi surfaces remain nondegenerate during 
the interpolation, we obtain the Z2 FSTI as the parity of 
the winding number given by Eq. 

N2D - (-l)"^"- = (-l)^2:.«gn(5.)Ci. ^ J] (isgn(5,))^^» . 

s 

Such a formula can be further simplified by noticing the 
following two properties of the Chern number Cis car- 
ried by the Fermi surfaces: i) The Chern number of each 
Fermi surface satisfies (— 1)*^!= = (—1)"*% where is 
the number of TRI points enclosed by the s-th Fermi 
surface, ii) The net Chern number of all Fermi surfaces 
vanishes, Cis = 0. We will leave a more detailed 
demonstration of these two conclusions to the supple- 
mentary material^, and only sketch the physical reasons 
for them here. The conclusion i) comes from the fact 
that a Fermi surface which only encloses one TRI point, 
such as Fermi surface 1 in Fig. [21 always enclose a sin- 
gularity at the TRI point due to Kramers' degeneracy. 
One can prove that the Chern number is always odd by 
making use of time-reversal symmetry. The Fermi sur- 
faces enclosing multiple TRI points can be adiabatically 
deformed into several Fermi surfaces, each enclosing a 
single TRI point. The conclusion ii) is a consequence 
of the Nielsen-Ninomiya theorem^S which states that the 
total chirality of a 3D lattice system must be zero. 

Using these properties of Cis, we finally obtain the fol- 
lowing expression for the Z2 FSTI which is independent 
of the interpolation to 3D: 

N2D = n MSs)r' ■ (9) 



4 



The criterion shown in Eq. ([9]) is quite simple: a 2D TRI 
superconductor is nontrivial (trivial) if there are an odd 
(even) number of Fermi surfaces each of which encloses 
one TRI point in the Brillouin zone and has negative 
pairing. 

Dimensional reduction to ID and generic ex- 
pression of the Z2 invariant. Following the same 
logic, the dimensional reduction can be carried out again 
to obtain the Zi FSTI in ID. This results in an identical 
formula to Eq. Since in ID each Fermi "surface" 

(which consists of two points at fc^ and —kp) always en- 
closes one TRI invariant point, the FSTI is simply 



ID 



(10) 



where s is summed over all the Fermi points between 
and TT. In other words, a ID TRI superconductor is non- 
trivial (trivial) if there are an odd number of Fermi points 
between and tt with negative pairing. Two examples 
wi □ 




FIG. 3: Simple examples of (a) nontrivial and (b) trivial pair- 
ing in a ID system. The red and blue dots are Fermi points 
with negative and positive pairing, respectively. The phase 
9 = d-nk is taken to be the same at fc = for the two bands. 
At fc = TT, ^ of the two bands are differ by 2-k for nontrivial 
pairing (a) and by for trivial pairing (b). 

Interestingly, from Fig. [3] we can get an alterna- 
tive understanding of the ID topological superconduc- 
tor, which can apply to a generic ID TRI superconduc- 
tor beyond the weak pairing limit. As discussed earlier 
in Fig. ([1]), the sign of the pairing Ss determines the 
winding of the phase 0„k across the Fermi point. On 
the other hand, we have shown that time-reversal sym- 
metry requires T^QuX — Q^]^, from which we can find 
that 9nk = 9n-k if \n,\i) and |n, — k) label a Kramers' 
pair. Thus, along the path from = to fc = tt, the 
change of Onk and 6n-k must be the same modulo 27r: 
Jp dk {dkOnk ~ dkOnk) = 27rri, n G Z. In the examples 



shown in Fig. [3] we have n = \ for the nontrivial pairing 
and n — Q for the trivial pairing. Such a parity difference 
of the winding number of 0nk turns out to be generic, 
and can be captured by the following Z2 FSTI: 



pf(rtQfc=o) 



exp 



dfcTr 



QldkQk (11) 



where we have used TtQ^T = Q'^^,^ ^ T'^Qk 

— {T'^Q-k)^ , so that T^Qk is anti-symmetric for k — 
0, TT, and the Pfafhan is well-defined. It is straightfor- 
ward to show that Nio = ±1 is a Z2 quantity, and also 
a topological invariant. More details on the properties of 
the Z2 FSTI dnil) and its relation to the FSTI ^ are 
given in the supplementary materials. Eq. (jDip is the 
topological superconductor analog of Kane and Mele's 
Z2 invariant in quantum spin Hall insulators^^ . Follow- 
ing the same approach as Ref.^^'^'', one can obtain three 
Z2 invariants in 2D, one of which is the "strong topolog- 
ical invariant" N2D = Nmiky = 0)NiD{ky = tt), with 
Nioiky = 0(7r)) the \D topological invariant defined for 
the ky = 0(7r) system, respectively. This topological in- 
variant is robust to disorder, and is equivalent to the one 
described by Eq. ([9]). 

In summary, we have presented the criteria for TRI 
topological superconductivity in the physical dimensions 
one, two and three. When the Fermi surfaces are non- 
degenerate, the criteria are very simple. In three dimen- 
sions, the winding number is an integer which is deter- 
mined by the sign of pairing order parameter and the 
Chern number of the Fermi surfaces. In one and two di- 
mensions, a pairing around the Fermi surface is nontrivial 
if there are an odd number of Fermi surfaces with a neg- 
ative pairing order parameter. We also obtained an ex- 
plicit formula for the Z2 invariants applicable to generic 
ID and 2D TRI superconductors. Our results provide 
simple and physical criteria that can be used in the search 
of topological superconductors. Our FSTI's suggest to 
search for topological superconductors in the nonconven- 
tional superconductors with strong inversion symmetry 
breaking and strong correlation. The strong inversion 
symmetry breaking is necessary to generate spin-split 
Fermi surfaces, and strong electron-electron Coulomb in- 
teractions prefer the pairing to have a nonuniform sign 
in the Brillouin zone/24 
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APPENDIX A: BACKGROUND AND DERIVATION OF THE BDG HAMILTONIAN 

We first list some basic properties of time-reversal invariant superconductors in generic dimensions and then go on 
to derive the form of Eq. (1) from the Letter. Consider a general TRI superconductor with the Hamiltonian 

H^Yl k'^k^k + I (^lAk^l^k + H.c.)] ^ Yl {^l ^"^k) Hik) ( % ) (Al) 



with 

"'^^ ~ ^ At(k) -/i^(-k) 
The normal state Hamiltonian ft.(k) is time- reversal invariant, which means there is a matrix T satisfying 

T-VkT = Tti/'-k, Tt/ikT = = -T^'^i T^^T = I (A3) 

From the transformation property of ^/'k we can obtain 

T- VjcT = ^l^T (A4) 

so that 

and the time-reversal symmetry of the Hamiltonian requires 

rtiJ(k)T = (-k)^, with T =('-'- j (A6) 



On the other hand, the following identity 

requires the particle-hole symmetry of the BdG Hamiltonian: 

Cti7(k)C = -H{--kf, with C = jj ^ j (A8) 

The two symmetries (|A6p and (|A8p require the pairing matrix A(k) to satisfy 

A(k) = -A^(-k), (rA^(k))^ = rAt(k) (A9) 



If we define 



X = *TCt = ( ) , (AlO) 



then we have 

xt-ff(k)x CTtiJ(k)TCt = Ci7(-k)'^C^ = -iJ(k) (All) 
The "chirality operator" x can be diagonalized by 

To derive Eq. 1 from the main letter we transform the basis to the eigenbasis of x to get the Hamiltonian form 



^(k) = Fi/(k)T/t - f _ ^^^t ''•^ ^ '^^^ ) (A13) 



k 
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As is mentioned in main text, the matrix = /ik+i^Aj^ can be decomposed by singular value decomposition (SVD): 

= h:, + iTAl = ulDi,Vi, (A14) 

in which Dk is a diagonal matrix with nonnegative real diagonal components, and f/k and Vk are unitary. The 
Hamiltonian -ff (k) can be diagonalized as 



H{k) 



t 



FjDkC/k 



^/2 V -^J jy -D^J V2\UI 



(A15) 



Thus we see that the eigenvalues of the Hamiltonian are given by the eigenvalues of and —Dk. For a gapped 

Hamiltonian, all the eigenvalues of Dk are positive, so that we can adiabatically deform to I, which deforms the 
Hamiltonian to the form 

#(k) -(^Qt^ ^ , Qk = UlVu e U(N) (A16) 

It can be seen from the derivation above that Qk is uniquely determined by the BdG Hamiltonian i?kj up to a 
k-independent U(N) x U(N) rotation 

Q^^gQkh, g,h€\J{N) (A17) 

All physical information carried by Qk, such as the topological invariants, is insensitive to this global U(N) x U(N) 
rotation. 



APPENDIX B: DETAILED DERIVATION OF THE 3D FERMI SURFACE FORMULA 

In this section we will show the detailed calculation of the 3d Fermi-surface formula. Beginning with the generic 
form of the winding number in 3d we will show how to derive Eq. 6 from the main letter. The general formula for 
the integer valued topological number is 



^w = ^J rf'ke^^^^TV [QidiQi^QldjQuQldkQu] ■ (Bl) 



First of all, if Ak = for some region of k, the winding number density vanishes in that region. To see that. 



notice that for Ak = 0, Qk is an adiabatic deformation of ^k = hu, so that Qu is Hermitian, and Q'l = Q^Qk = I- 



Consequently the winding number density is given by 

1 



Pw 
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e^^'^Tr [Qk^^QuQl.^JQkQl.^kQu] ■ (B2) 



By making use of di (Qk) — QkdiQk + diQ^Qk = 0, i.e., {Qki^^Qk} = 0, one can prove that pw = 0. This 
confirms our statement that in the weak pairing limit, when only the pairing around Fermi surfaces is considered, the 
topological invariant Nw is completely determined by the physics in the neighborhood of the Fermi surfaces. 
As discussed in Eq. (4) of the letter, in the weak pairing limit Qk can be written as 

Qk = 5Ze'^"''|n,k)(n,k| (B3) 

n 

with e*^"'' = (e„k + i^nk.) / l^nk + i^nk\ and dnk = {'n,'k\ TA]^ |«, k). To the leading order, near the Fermi surface we 
have 
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In the limit (5. 



0, we have 



hm 6'„k —f vrsgn {5nkp)ri {kp - k±) 

with r]{x) the step function satisfying ri{x) — I, x > and 77(2;) =0, x < 0. Thus we obtain 

dk^Oni^ = -7rsgn((5„k)'5(A:j^ - fcf) 
In the vector form, this equation can be written as Eq. (5) of the letter: 

V6'„k = -7rv„kSgn ((5„k) S (e„k) 
Now we simplify the winding number formula by using Eq. (|B3p . After some algebra we obtain 



(B5) 



(B6) 



(B7) 



Nw = 



27r2 



2i 



9,0„(af sin^ 



at sm — — 



sp ■ "sp 



(B8) 



where 9ns = 0n — (^s and a"'* = —i {n, k| 9^ |s, k) is the non-Abelian adiabatic connection. When we restrict the pairing 
to an energy shell — e < e„k < e and take the e — > limit, the only nonvanishing term is the one with di9„ which has 
a 5 function on the Fermi surface. This leads to 



w 



'2^ 



r pkp+e/vp n 

/ d^kii / dfc_L V9i6'„sin2^(arar-a2'ar) 



' kp—e/vp 



dda sin^ — 



FS 



( ns sn ns sn\ 



(B9) 



where 9^ are the values of 0„ right outside and inside the Fermi surface, respectively. When there is only one band 
that crosses the Fermi surface, 9^ = 9^ for all other bands. Labelling the single band crossing the Fermi surface with 
n = 0, we have 



Nw 



'2^ 



FS 



% - sin (00 - Os) 



(BIO) 



Since 9^ and 9s all have the values or tt, the second term sin(6'o — 9s) vanishes, and we have 



w 



1^ 



' s^O 

= ^^sgn(<53k) / d^k||(9ia°°-a2a;'') 

^ FS ''^^ 

= ^E^Sn('^s)Ci^. 



(Bll) 



It should be noted that, (i) The superconducting gap on the Fermi surface is given by |(5„k|, so that sgn((5„k) is the 
same for all the k on the same Fermi surface, otherwise the superconducting gap would vanish for some k. (ii) The 
Chern number of the s-th Fermi surface Cis is defined with the normal vector along the direction of v^, which is 
opposite for an electron pocket and a hole pocket. 
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APPENDIX C: PROOF OF CHERN-NUMBER PROPERTIES USED IN THE DIMENSIONAL 
REDUCTION TO 2D SECTION OF MAIN LETTER 

In the main letter, we have used the following two properties of the Fermi surface Chern number to obtain the 2d 
Z2 formula: 

1. The Chern number of each Fermi surface satisfies {—Vf^^' = (—1)™% where nig is the number of TRI points 
enclosed by the s-th Fermi surface. 

2. The net Chern number of all Fermi surfaces vanishes, C\s — 0. 
In this section, we will prove both properties. 

1. Proof of Property 1 

We first study a simple Fermi surface enclosing one TRI point, e.f?., the F point, as shown in Fig. 2] (a). Denote 
the states at the Fermi level as |s, k). The Berry phase gauge potential is defined by = —i (s, k| di |s, k). In the 
following, we will denote ai = af^ for simplicity. The time-reversal invariance of the normal state Hamiltonian 
requires the time-reversed state T (|s, k)) to also be on the Fermi surface. When the bands are non-degenerate on the 
Fermi surface, in general we have 

r(|s,k» = e^^-\s,~k) (CI) 

^a,(-k) = -»(s,-k|^^|s,-k)=iT((s,k|)e''^''9,[e-''^''T(|s,k))] 

= di(fi^ + i{{s,k\di{\s,k)))* 

= d^ipk + a,(k). (C2) 

Thus the gauge curvature is 

/y(k) = a,a,(k) - dM^) = (C3) 
We denote the upper half of the Fermi surface with A:^ > as FS+ and the lower half as FS_. Thus 

Jfs ^'^ Jfs+ Jps_ 

Since the two form dfl^^ denoting the normal direction of the Fermi surface is also odd in k, the contributions of FS+ 
and FS_ to the Chern number are equal, so that 

Cis = - I rf^^'V^J(k). (C5) 

Jfs+ 

Since FS_|_ is a manifold with boundary, the Chern form is equivalent to a boundary integral: 

Cis^- (f dl'aiik) (C6) 

JdFS + 

in which 9FS+ is the boundary of FS+, i.e., the = section of the Fermi surface, and dP is the tangent vector to 
(9FS+. However, it should be noted that Eq. (jC6p holds only if ai{k) is continuous in the whole FS+. In a generic 
gauge transformation ai(k) — > ai{k) + di(f)\^ on the boundary 9FS+, the right hand side of Eq. (|C6p can change by 
an even number: 

- (f dl'ai{k) -^-i dUa^{k) + - (f drd,(j}]^ = - / dVai{k) + 2n, n e N (C7) 



Thus we have 



Cis^- I dPai{k) mod 2 (C8) 
JaFS+ 
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in a generic gauge choice. 

Since the section fc^ = of the Fermi surface is also symmetric under time-reversal, we can split it to two parts Li 
and L2, which are the time reverse of each other, as shown in Fig. |4] (a). Noticing that the tangential vector d/* is 
opposite for k and — k, and by making use of Eq. (|C2p we have 

® drai(k) = - ® dl^aiCk) — (p dPdi(p\^ 

^ Ci, = - i dPa,(k) = -( (f + / ) dPa,{k) =-- (f dPd^^k. (C9) 
I" JaFS+ \Jl, Jl^J ^ Jl2 

One can always split the boundary so that there are only two points A and B on the interface between Li and 
Due to time-reversal symmetry, the two points must be the time-reversed partners of each other, and the formula 
above becomes 

Cis^--{Va-^b) mod 2 (CIO) 

TT 



Denote the momentum of A and B as k^ and k^ = — k^i, according to the definition Eq. jCTj) we have 

T{\s,kA)) = e*'^-^|s,kB), r(|s,kB)) -e*'^- |s,kA) 
^T{T{\s,kA))) = T(e*'^-^|s,kB)) =e-"^-^e*'^-|s,kA) (Cll) 

On the other hand, we have = — 1 for each state, so that 

„i(VA—VB) 1 ^ 1 ™„J o (C12) 

Thus we have proved t 




FIG. 4: (a) Schematic picture of a Fermi surface enclosing one TRI point (0, 0, 0). The Fermi surface is separated to two parts 
FS+ and FS- hy — plane. The interface between the two parts is further split into curves Li (red curve) and L2 (blue 
curve) which are time-reversal partner of each other. The interface of Li and L2 are given by points A and B. (b) Schematic 
picture of a Fermi surface enclosing two TRI points (0, 0, 0) and (0, 0, vr). Similar to (a), the Fermi surface is separated to FS+ 
and FS-, and the interface between FS+ and FS- are split into Li (red curve) and L2 (blue curve), which intersect at two 
pairs of points Ai,Bi and A2, -B2. 

For the Fermi surfaces enclosing more TRI points, as shown in Fig. S] (b), the proof is similar. Due to the time- 
reversal symmetry, we can always reduce the Chern number to an integral over the upper half of the Fermi surface 
FS+ as in Eq. (jC5|) and (|C6|) . Generically, FS+ has two boundaries at fc^ = and — tt, so that 



TT 




Cis = ^ f dn'^f.j (k) - - I i - f I dPa^ik) (C13) 

'fS+ ^ \Jd„FS+ JdoFS-i 



where 9o,7rFS+ stands for the boundary of FS+ at fc^ = and k^ — n respectively. In the same way as above, 
the boundary at fc^ = can each be separated into two parts Li and L2, with several pairs of interface points 
Ai, Bi, i = 1, 2, ..,po- By the same derivation as above one can prove ipAi — = mod 27r, and 

If 1 

-d) dPa,Ck) = V((^A. -m) =Po niod 2 (C14) 
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The same argument works for the fc^ = vr boundary. Denoting the number of interface points at the fc^ = tt boundary 
by , we have 

Cis =P7v -Po mod 2. (C15) 

If the boundary 9o,7rFS+ encloses ttio^tt number of TRI points, respectively, we have pq = rriQ mod 2, p^r = 'mTr mod 2. 
Thus 

(_l)Ci. ^ (_;l)P,-P0 ^ (_;L)m,+mo ^ (C16) 

with rUs = + mp the total number of TRI points enclosed by the Fermi surface. Thus we have proved the property 
1. 



2. Proof of Property 2 



To prove property 2, we take a simple s-wave pairing 

Ak = AqT (C17) 

with Aq a real number. For such a pairing the matrix ylk = /ik + ^'^Aj^. = /ik + iAgl, so that the pairing on all the 
Fermi surfaces has the same sign: 

<5, = (s,k|rA[|.s,k) = Ao, Vs. (C18) 
Consequently, the topological invariant in the weak pairing limit is given by 

N^{Ao)^^J2'gniSs)Cu = '^^^^J2^is. (C19) 

s s 

On the other hand, the BdG hamiltonian (IA13I1 for this simple pairing can be diagonalized easily to obtain the 
eigenvalues 



^„k = ±V<k + A^. (C20) 

Thus for finite Aq, the spectrum of the BdG Hamiltonian is always gapped, so that the winding number Nwi^o) 

remains invariant for all Ag > 0. Thus we can compute Nw in the limit Aq +00. The unitary matrix Qk is given 
by 

^ lim Qk = iS^ \n,k) {n,k\ ^ il (C21) 

n 

Obviously, the winding number Nwi^o ~^ +00) — 0, so that Nwi^o) — for any Aq. According to Eq. (|C19[) we 
have proven property 2: 

Cis = 0. (C22) 



APPENDIX D: PROPERTIES OF THE ID Z2 TOPOLOGICAL INVARIANT (11) 



In this section, we will study some basic properties of the Z2 topological invariant defined in Eq. (11) of the letter, 
and show how it is reduced to the Fermi surface formula (10) in the weak pairing limit. 
We start from Eq. (10) of the letter: 



Id 



Pf (TtQfe^o) 



exp 



dfcTr 



QidkQk 



(Dl) 
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First of all, T^Qk is antisymmetric since 



rtftwT = h^^, rt ( TAl) T = At r 



TAT 



(D2) 



Thus the Pfafhan is well-defined at fc = and k — it. 

Since Qfc G U(N), we have det Qfc — e'^*' which is a U(l) phase. Since Tr Q\dkQk = Tr[logQfc] = logdetQ^ — 



iipk, we have dkTr 



QldkQk 



exp 



i {(p{tt) — (p{Q)) mod 2tt, so that 



dkTv 



QldkQk 



Thus 



det (rtQfc^,) 
= det(rtQ,.o) 



dfcTr 



dot (rtg^^o 

det (TtQfc^, 
Q^a^Q^. I = 1 



so that Nid always takes the value of ±1. 

Now we show that iVi^ is a topological invariant. For an infinitesimal deformation Qk > Q'f. = Qk 



factor exp [-^ dkTv QldkQk 

dkTi 



only depends on the deformation of Qk at fc = and n: 



exp 



Q'kdkQ'k 



exp 



dfcTr 



QldkQk 



e 2 



On the other hand, the change of Pfafhan is given by 

Pf {T^Q'k^o,.) = ei^'^-°-Pf (rtQUo, J 



(D3) 
(D4) 

5Qk, the phase 
(D5) 
(D6) 



Consequently, we see that 6Nid = in any smooth deformation of the unitary matrix Qk as long as time-reversal 
symmetry is preserved. 

In the weak pairing limit, the general formula (jPip can be reduced to the Fermi surface formula given by Eq. (10) 
of the letter. In the weak pairing limit, assume there are M Fermi points kg, s = 1, 2, AI between and tt. As 
discussed in the main text, we require that the Fermi level does not cross any band at /c = or tt. As discussed in 
Fig. 1 of the letter, each Fermi point leads to a domain wall of Ogk for the corresponding band s crossing the Fermi 
level. According to Eq. (|B3[) in the weak pairing limit we have 



det 



Qk = exp i ^ Onk 

\ n 



(D7) 



Across each Fermi point kps, the phase 6sk will jump by — Trsgn (vps^sk^ ) and the 9nk for other bands remain invariant. 
It should be noted that the sign of vp enters the expression since the winding of 9sk is given by — Trsgn (J^fc^) along 
the direction of the Fermi velocity vpg. Consequently, the phase log det Qk = i ^"-^ is changed by — iTrsgn (vpsSsk^) 
across the s-th Fermi point, and the net change of log det Qk from to tt is given by 



exp 



„7r HI 

/ dkdk log det Qk = -iir ^ sgn ( 



dfcTr 



QldkQk 



VpSskJ 

S S 

Yl (sgn (SskJ) Yl (-isgn (vps)) 



(D8) 



When there are m Fermi points with positive vps and n Fermi points with negative vps, n — m gives the number of 
bands which are above the Fermi level at fc = 0, but below the Fermi level at fc = tt. If we denote A''2(0) and N2{tt) 
as the number of bands occupied at fc = and k — n, respectively, then n — m ^ N2{tt) — N2{0). Since all bands are 
paired in Kramers pairs at fc = 0, tt, A'^2(0) and N2{t^) must be even. Thus we have 



62 



7r{n— m) 



(-1) 



W2(^)"W2(0) 



(D9) 
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Now we study the Pfaffian Pf (T^(3fe=o.7r) • Since we have assumed the Fermi level does not cross the bands at 
/c = 0, TT, in the weak pairing limit we have A/j^o,7r — 0- If the normal state Hamiltonian is diagonalized to 



hk = Ul 



eAr(fc) 



(DIO) 



Qk can be obtained by 



-I 



N2XN2 



u. 



(Dll) 



in which iVi and N2 are the number of unoccupied and occupied bands, respectively. The Pfaffian Pf (T^Qk) can be 
obtained as 



-I 



N2XN2 



Uk 



Pf ( (jfrtt/M ^^^x^i 



XN2 



■ det Uk 



By making use of the time-reversal invariance condition T^Q^T = Q^^, one can prove that 



IIJV2 X N2 



= 



(D12) 



(D13) 



for fc = 0, TT, so that the matrix ul'^T^U'l is block diagonal. Consequently, we have 



Pf ( f/frtc/t 



mi X Ni 



IWoxAfo 



(_l)JV./2pf |^[^tT^tt/t) = (-l)^^/2pf (rt) .dctt/2 



(D14) 



Thus 



pf(rtgfe) = (-i)^^/2pf(rt) 



(D15) 



for k — 0,Tr. It should be noted that the number of occupied bands N2 is always even for k — 0,7t due to Kramers 
degeneracy. 

Combining Eq. (|D8|) . (|D9|) and (|D15p we obtain 



(DI6) 



Thus we have proved that the general Z2 invariant (11) in the main text is equivalent to Eq. (10) in the main text in 
the weak pairing limit. 
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